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laws in the classical Toda ring. We also show that perturbing the integrability leads to a decay of 
the currents at long times, with a time scale that is determined by the perturbing parameter. We 



1-^ ' summarize several known results concerning the Toda ring in 1-dimension, and present new results 
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relating to the frequency, average kinetic and potential energy, and mean square displacement 
in the cnoidal waves, as functions of the wave vector and a parameter that determines the non 



I. INTRODUCTION 

Toda's non linear lattice [l| is one of the very few examples of non linear lattices in 
condensed matter physics, where explicit analytical solutions are available for the dynamics. 
There are several aspects of condensed matter physics where the Toda lattice is a useful 
model. Toda himself applied his non linear lattice to understand heat propagation 2|, 
and further studies with added impurities throw interesting light on this phenomenon 3|. 
Interestingly, the lattice has also found recent applications in the context of the dynamics of 
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DNA j^, where the Toda interaction is a reasonable representation of the known non-linear 
couplings between base pairs. It has also been used to represent the potential of hydrogen 
and peptide bonds in the a- helix |5[. 

In this paper we study in detail the Toda lattice with periodic boundary conditions (the 
Toda ring). Our aim is two- fold. Firstly we derive several new results. Secondly, since the 
Toda lattice and its properties are less well known to students of condensed matter than 
they deserve to be, we collect together some of the basics of the model and its solution in a 
form and notation that is standard in condensed matter physics. 

The excitations of this lattice are not phonons as in a harmonic lattice, but can be 
expressed in terms of non linear excitations that are termed solitonic. For periodic boundary 
conditions the excitations are more properly the cnoidal waves corresponding to a family of 
waves characterized by a wave vector and another parameter, related to the non linearity of 
the excitations, namely the elliptic parameter m discussed below in Sec. IIIAI The lattice is 
very simple to describe and the solution is both beautiful and instructive. Surprisingly, we 
find that the dispersion relation of the excitations given by Toda is only correct in the limit 
of weak anharmonicity or long wavelength. Here, we give a complete expression, which does 
not seem to have been calculated before. 

Soon after Toda found the exact solution, his model was found to be exactly integrable [g, 
|7|, i.e. it has an infinite set of "generalized conservation laws". These conservation laws are 
expressed through conserved currents that Poisson commute with the Hamiltonian as well 
as each other, and the stability of the solitons is understood to arise from the existence of 
these currents. 

There is considerable interest in the role of the conservation laws in the transport of heat 
or energy, and the Toda lattice provides an excellent model to test some ideas about their 



role in transport, as detailed below. Since the model is classical we can study reasonably 
large systems (up to 64 atoms in our largest studies), unlike in quantum systems where the 
Hilbert space grows exponentially with the number of sites. 

We address two specific issues in transport theory in this paper. The first is the role 
of integrability of the model in determining the exact value of the asymptotic correlation 
function of the energy current — this value provides us with the coefficient of the delta func- 
tion in the thermal conductivity at zero frequency. It is known as the Drude term in the 
Kubo conductivity, and is widely discussed in current literature. The second issue concerns 
the role of perturbations of integrable models, whereby conservation laws are destroyed. 
We present results on the decay of the energy current in a slightly perturbed Toda lattice 
for various values of the parameter that destroys integrability, and show that there is an 
underlying scaling picture which provides a general understanding of this phenomenon. 

The plan of the paper is as follows. In Sec. |TT]A we define the Toda ring, and in Sec. HTlB 
we discuss the extreme limits of harmonic and anharmonic interactions. The frequency of 
periodic solutions (cnoidal waves) of the Toda ring are derived in Sec. HIlC. In Sec. HTlD we 
highlight the differences between the results for the frequency spectrum derived in Sec. Hi] 
C and the spectrum determined by Toda. We also show how, in the extreme anharmonic 
limit and at long wavelength, the cnoidal waves can be viewed as a train of isolated solitons. 
Next, in Sec. Illlt we calculate the kinetic energy in the modes as a function of wave vector 
and anharmonicity parameter, and also the ratio of the kinetic energy and the average 
displacement as functions of wave vector and anharmonicity. In Sec. [IV]we list the conserved 
currents of the Toda ring obtained from the Lax matrix. Next, in Sec. |Vl we analyze the 
persistence of energy currents by expanding the energy current in terms of these conserved 
currents. Comparing with numerical results we show that the persistent part of the conserved 
energy cannot be expressed in terms of the Lax currents alone, but quadratic combinations 
of Lax currents are also needed to get an accurate description. In Sec. IVll we consider 
a Toda ring in which a small interaction is added which breaks integrability. We calculate 
numerically the decay of the persistent currents as functions of time for different values of the 
perturbing parameter. We show that the results fit a scaling property that is expected, but 
with substantial corrections to scaling so very long runs are needed to reach the asymptotic 
scaling regime. Finally, in Sec. IVIIt we summarize the main results of the paper and comment 
on them. Appendix |X] contains some needed results on elliptic functions. In Appendix [B] 



we discuss Toda's result for the frequency spectrum wj of periodic waves, and explain why 
this result is only correct in some limiting cases. Appendix [Cl gives an alternative derivation 
of the dispersion relation of the cnoidal waves, while Appendix [D] computes the potential 
energy of the Toda ring. 

As discussed above, one of the objectives of this paper is to collect results on the Toda 
lattice in a form accessible to condensed matter physicists. In addition, our result for 
the dispersion of cnoidal waves in Eq. ( !T6|) has not been published before, and the fact that 
Eq. (TT6|) differs from the result given by Toda, is new. Other new results are the computation 
of the average potential and kinetic energies in the cnoidal waves in Section IIIH and the 
expansion in Eq. fl23p . Finally, our results for the connection between persistent energy 
currents and the Lax conserved currents in Sec. El and for the decay of the perturbed Toda 
ring in Sec. IVIt are new. 

II. TODA RING, CNOIDAL WAVES AND THEIR SPECTRUM 

We summarize some interesting facts about the Toda lattice in this section. By imposing 
periodic boundary conditions, we deal with a ring of finite extent. The Hamiltonian of the 
Toda lattice is 
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The displacement variable m„ is defined through i?„ = i?° +-u„, where i?° is the equilibrium 
position of the n^^ atom. In principle the displacement u„ ranges between ±oo, so for 
consistency, we must either imagine that the lattice constant -R^+i — -R° is infinite as well, 
or that the displacements are transverse to the ring. Further we assume m„ = Un+N as 
appropriate for a ring geometry. The two non-exponential terms in the interaction potential 
are irrelevant, since they add but a constant to the energy, but it is convenient to include 
them since the "two body potential" then explicitly displays a minimum at zero relative 
displacement. The variable p„ is conjugate to Un satisfying the standard Poisson Bracket 
(PB) relation {un,Pm} = Sn,m, and a,b,M are parameters. 
The equations of motion follow from Hamilton's equations. 



We note that the total momentum ptotai = ^nPn is a constant of motion. The dynamics 
has in fact many more conservation laws, a consequence of the property of integrability 
which was proved by Henon |6| and Flaschka [3] for this system. The explicit form of 
the conservation laws are given later in the paper, and follow from the Lax structure that 
underlies the dynamics. The parameters a, b, and M give us explicit freedom to interpolate 
between the harmonic and extremely anharmonic limits. 

A. Harmonic and Anharmonic Limits 

A formal Taylor expansion of the exponential interaction gives us 



No N 



n=l n=l 

where k = ah. As long as the displacements satisfy \uj\b ^ 1, the anharmonic terms do 
not become important, so we expect a harmonic response. However for \uj\b 3> 1, the 
anharmonic terms will dominate. To recover the harmonic lattice, one could formally take a 
limit 6 — > and simultaneously let a — )■ /t/fe with k remaining finite, so that the anharmonic 
terms are explicitly killed. 

We will see below that the cnoidal waves of Toda, contain a parameter, the elliptic "m" 
parameter, which controls the amplitude of the waves, and varying this parameter gives 
harmonic as well as strongly anharmonic response. Physically, the elliptic parameter "m" 
may be viewed as tuning the anharmonicity, with the harmonic limit being m — )• and 
the extreme anharmonic limit being m — )■ 1. Mathematically, the parameter m plays a 



fundamental role in the theory of Jacobian elliptic functions[l8|. 

For later use, we note that the harmonic limit has a dispersion ujk and sound velocity cq 
given by 

k 

Uk = Co 2|sin-|, 




CO = \TJ- (4) 



B. Single parameter formulation 

The Toda lattice is also integrable in quantum theory. Quantum integrability was estab- 
lished by Sutherland |9|, Gutzwiller 10[, Sklyanin |ll|, and Pasquier and Gaudin 12[ using 
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who also establish their equivalence. In the viewpoint of 

lattice emerges from a crystallization of a gas of impenetrable particles with an interaction 

oc {sinh(i?„-i?„)}-2. 

Quantum mechanically, it is possible to reduce the Toda lattice to a single parameter 
problem [8J. We scale m„ — )■ 6 «„, Pn ^ \ Pn, H — )■ jjjjH and set rj = 2Ma, so that 

N 2 ^ 
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In this representation r^ — )■ oo gives the harmonic limit, since the displacements become 
very small so the potential energy remains small. On the other hand, ?7 — )■ corresponds 
to the extreme anharmonic limit, since now the displacements are large, and hence high 
order terms in the expansion of the exponential matter, and ultimately dominate. We need 
to keep in mind that in this extreme non-linear limit of the model, the particles are not 
allowed to cross so they act as impenetrable billiard balls |8| with free propagation between 
successive collisions. Classically, even this one parameter rj can be removed by a rescaling 
of the displacements, so large anharmonicity corresponds to large displacements and vice 
versa. 

C. Cnoidal wave solutions 



We derive here the formulas for the excitation spectrum of the Toda ring. The relation 

His papers l|, |2| and book 13[ focus on 



with Toda's work is discussed in the Appendix 
a set of dual variables and give a solution for the displacement, but we point out that his 
dispersion relation is not appropriate for the periodic boundary conditions (i.e. a ring) which 
is the focus of this paper. Here we present an explicit solution for this case, which is not 
available in literature as far as we can tell iJ] . Our Eq. ( HTll below, concerning the total 
energy to mean square amplitude ratio, and related results, also seem to be new. We feel 
that they, too, are helpful in appreciating the Toda system from a condensed matter point 
of view. 

The Toda Hamiltonian in Eq. ([1]) leads to the following equation of motion for the 
displacement J28|| : 
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We seek a special kind of solution namely a constant profile solution 

b Uj{t) = dk{<pj{t)), 
4>j{t) = k j -ujkt = k{j - Ckt) = 27r ( - - z/t j , (7) 

where (pjit) is the usual phase factor depending linearly on space and time with wave vector 
k (or equivalently wave length A), and angular frequency Uk (or equivalently frequency z/), 
and we have defined the velocity by Ck = uJk/\k\. We will often omit the argument of the 
phase for brevity. We set 

k = I-, (8) 

where the A^ integers u obey 1 < u < N. Consequently, the Brillouin Zone is the range 
< A; < 27r. Periodic boundary conditions, m„ = Un+N, are satisfied if the function dk is 
periodic, i.e. dk{x) = dk{x + 27r). Here Uk is yet to be determined, along with the form 
of dk- We now define a scaled frequency ujk by Uk = \/jj ojk so the equation of motion, 
Eq. ([6]), reduces to a non-linear, differential, difference equation; 

col dl{<p) = e'^*(^-'=)-'^^(^) - e^'=W-'='*(^+^). (9) 

This equation is satisfied by the choice 

where we summarize, in Appendix [Aj the necessary definitions of the elliptic theta functions 



as needed for this work |29|. To see that Eq. ( ITOj) solves Eq. (J9]) we use the Jacobi addition 
formula for the 9 functions [l5|, and write the first term in the RHS of Eq. (jQ]) as 

^dj^(<l>~k)-dk((t>) _ 




(11) 

The second term in the RHS of Eq. ^ is obtained by replacing by + /c. Upon using the 
relationship between the theta functions and the Jacobian elliptic functions [16| 



OAx) 1 ,2K , 

-— - = m4 sn X , (12) 

6Ax) TT 



we find the RHS of Eq. (Q is given by 
RHS Eq. dHD 
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where K = K{m) is defined in Eq. (jM]). The LHS of Eq. (E]) is 
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Thus Eq. iQ is satisfied provided we assign the frequency as one of two equivalent expres- 
sions, 
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with the help of the following standard relations 17[ among Jacobi's constants: 



^'i(O) = ^2(0)^3(0)^4(0), ^2(0)^3(0) = -m^ K. 
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Equation fITBl) does not appear to have been published before, but an equivalent result has 
been obtained by Sutherland [ij] in unpublished notes. 

To summarize this section, the displacement given by the expression in Eq. (fTOl) satisfies 
the equation of motion, Eq. iQ, for the Toda ring with the frequency given by Eq. flT6|) . 



D. Connection between the frequency spectrum of the cnoidal waves and Toda's 
spectrum. 



Toda's result for the frequency of wave-like solutions of the Toda lattice is 
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We give a derivation of this in Appendix |Bl The difference between our result for the 
dispersion relation in Eq. flT6|) and Toda's in Eq. flTSj) is shown in Fig. [H For small m these 
are very close indeed, whereas for m — )■ 1 our solution always has a higher frequency. The 
difference between these dispersions also plays a role in determining the average potential 
energy as we show below in Eq. f l37p . As pointed out in Appendix [Bl the dispersion wj 
comes from a calculation which does not respect the boundary conditions, and hence the 
discrepancy with our result is not surprising. Perhaps what is surprising is that they are so 
close for small k, and differ significantly only for m close to unity. 
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FIG. 1: The main figure shows the dispersion relations at tti = 0.9999, which is strongly anhar- 
monic. The inset is for m = 0.7, which is only moderately anharmonic. The upper curve in both 
cases is the true spectrum, ZJ^, from Eq. (fT6l) and the lower one is tj|' from Eq. P^ . Except for m 
very close to unity, the two spectra are very close. 



In order to appreciate the role of the non linearity, we may usefully express the dis- 
placement and spectrum as expansions in power of the elliptic nome parameter q defined in 
Eqs. f lATllAal) : 



|2sin(|)| [1 + Vsin2(|) ^ Ug^ .i^^^^)] + o{q'), 



ul = |2sin(-)| [1 + Vsin2(-) cos(A;) + gSin2(-){6 - 5cos(A;) + 14cos(2A;) - 3cos(3A;)}] + O(g^), 



dk{<P) = 2q[cos{(f)) — cos{k — (p)] + 2q'^ sm{k) sm{k — 2(f)) 

+^q^ [cos^(0) - cos^(A; - 0)] + g^ sin(2A;) sin[2(A; - 20)] + O(g^). 



Here we made use of the expansion of 1/sn^ in example 57, page 535 of Ref. 18|] to rewrite 
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Eq. (fT8|) . Note that the difference between the first two hnes appears at 0{q'^), and is 
amphfied for large k. 
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FIG. 2: The displacements at a fixed time for different sites obtained from Eq. (jlOp . The data is 
for m = 0.9999, which is strongly anharmonic, and m = 0.7 where the anharmonicity is weaker. 
Here b = 1 and the ring is of length A^ = 100. In each case, results are shown for the two smallest 
values of k, see Eq. ([5|). 



In Fig. (|2]) we plot the displacements at different sites, obtained from Eq. (fTOj) . for a ring 
of length A^ = 100 at two different values of the parameter m with the lowest two A;-values, 
see Eq. ([8]). Notice the asymmetric shape of the displacements and their relatively broad 
structure. 

To further visualize the periodic solution, we first note that Eq. (ITOl) is closely related 
to the singly periodic Jacobi zeta function Z(u) = Z{u + 2K) via the relation ^dk{<p) = 
— [Z( — (0 — k)) — Z(—(f))]. The Jacobi zeta function has a formal Fourier series expansion, 
as well as one in terms of tanh. These are 
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tanh{^(«-2irz/)}, 
where the first (second) series expansion is particularly useful in the limit m ^ {m ^ 
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1) 27|. The solution of the Toda lattice on a ring, i.e. the "cnoidal wave", corresponds to 
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The last line follows from using the Poisson summation formula to rewrite the second line 27 1. 
Let us comment on the limit when the elliptic parameter m tends to 0, where the second 
line in Eq. f l23|) is useful. It is easily seen that we obtain the harmonic excitations in this limit. 
Using the standard expansion of the various objects in m, the displacement is given by u^ = 
^(cos(0j) — cos(0j — k)), and the prefactor of m makes these small amplitude oscillations. 
The spectrum is also that of the harmonic limit, since Eq. flT6|) becomes w^ = 2| sin|| and 



the phase factor in Eq. ([7]) is given by 0j = k{j — Ckt) with velocity c^ = ^Ja h/M = cq, as 
expected from Eq. ^. 

In the other extreme limit of the elliptic parameter, m — )■ 1, the elliptic functions degen- 
erate into hyperbolic functions and the cnoidal wave is regarded as a train of solitons, so 
that the periodicity of the displacements around the ring is unimportant. This is illustrated 
in Fig. (|3]), where we plot, for a typical case, the displacements Uj, their nearest neighbor 
differences Uj — Uj+i, and the exponential of the latter, i.e. exp {uj — Mj+i}. 

Let us now extract a single isolated soliton from the solution. We start with Eq. ( IDip 
for the exponential of the displacement difference in terms of sn^ and using dn [u) = 1 — 
m sn^(M), Eq. (1221) and the relation -^Z{u) = dn^(u) — E/K, we find a useful and formally 
exact series representation [jj, |2|, |l3| 

K °^ K 

^ + ^1 (^)' E sech2«-(0.(i) - 2^^), (24) 
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where a = ^ and 
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The periodicity in the phase angle — )■ + 27rz/ is manifest in this way of writing the 
displacement difference. 

Let us focus on 771 — )■ 1 so that we may set K' — )■ 7r/2 and a — )■ 1. With (f)j{t) = kj — 
Cq cJfct, we observe that the separation Aj between peaks of the exponential of displacement 
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FIG. 3: The isolation of a soliton from the cnoidal wave is ihustrated here. We plot various objects 
for m = 0.9999, 6 = 1 and a ring length N = 100 with k = Att/N . In the left panel we display the 
cnoidal wave displacements Uj, and the difference between successive displacements Uj — Uj+i. In 
the right panel, we display the exponential of the displacement difference. The latter shows two 
Toda solitons, with a clear separation ~ N/2 and a width that is several lattice constants, but 
much less than the separation. Note that while the displacements appear delocalized and wave 
like, the solitons are quite localized and hence particle like. 



difference Eq. (12^ . is 2t: /k. We shall see in Eq. (!32|) below that the width of these peaks is 
given by K^^ where 

(25) 



K = — k K, 



and so the requirement that the separation between the peaks is much greater their separa- 
tion is 

K ^ 1, or equivalently 1 — m ^ 1 . (26) 

It is clearly necessary that the width of the peak is at least several lattice spacings, and so 
we also need the condition k ^ 1. Combining with Eqs. (l25l) and ( l26l) the condition for the 
oscillations to be described by well separated solitons is 



k<K-^<€. 1, 



(27) 



which, as shown in Eq. ( !26l) . implies that m is very close to 1, i.e. the system is in the 
extreme anharmonic limit. 



In the proximity of a peak, we drop the sum over z/, set A — )■ 1 Si], and write 
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From Eq. (fT6|) . 6i{k/2) determines the dispersion, and in the expression Eq. (lASp . the first 
few terms in an expansion for positive k read as 
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As m ^ 1, K is large and K' — )■ 7r/2, and in the regime kK ^ 1 which is well satisfied here, 
see Eq. ( !27ll . we may further approximate this by writing 

71 , fkK\ , , 

a;fc = - smh (^— J . (30) 

The same answer is also found from Toda's relation Eq. ( ITSj) . since k is small enough in 
this regime that the distinction between the two dispersions is negligible. In terms of the 
parameter k defined in Eq. ( I25l) above, and calling the soliton velocity 

Ck = Co , (31j 

where Cq = a/||, we find 

gbK-n,+0 _ ^ ^ pl± _ (32) 

cosh [k(j - cj)\ 

Equation (|32|) is the profile of the famous Toda soliton, two of which are seen in the right 
panel of Fig. O We see that the soliton parameter k controls the amplitude of the solitons 
(through the prefactor in the amplitudes sinh^(/T;)), the velocity c^ and also the length scale 
of spatial variation ^. Since ^^^ > 1, the velocity of the soliton c^ in Eq. ( 1311) is always 
greater than the sound velocity Cq. 

For completeness we note the displacement for an isolated soliton, by using the same 
limits as above in Eqs. (126||27|) . The displacement Uj and its time derivative can be written 
from Eq. ( l23l) . where we retain only z/ = (since K ^ 1) and write 



Uj ~ - log 



cosh k{j — 1 — Ckt) 
cosh k{j — Ckt) 



(33) 



iij ~ [tanh k{j — Ckt) — tanh k{j — 1 — c^t)] . (34) 

It is now straightforward to see that the equation of motion Eq. ([H]) is satisfied by the single 
soliton solution given in Eqs. ( I32l) -(1M1). Toda notes that if we ignore the periodic boundary 
conditions, and imagine Eqs. (!32|) -( l34|) to be extended for all — oo < j < C)0, then integrating 
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Eq. (l33l) gives Mqo — U-oo = X' ^° there is a net compression near the sohton. Hence the 
sohton can be regarded as a local compression propagating with speed c^ given by Eq. fl3T|) . 

The momentum Miij vanishes at j — )■ ±00, as does the potential energy term Eq. (1321) . 
Therefore the energy of a soliton, obtained by ignoring periodic boundary conditions and 
opening the ring into an infinite chain, is finite. Substituting for iij and the displacement 
difference from Eqs. (133|) -( 1M|) into the Hamiltonian Eq. ([1]), and with pj = Miij, we obtain 
this energy to be esouton = ^(sinh Kcosh n — k). 

With the Fourier expansion in Eq. ( 123|) one can give an alternate derivation of the dis- 
persion relation of the Toda lattice with periodic boundary conditions in Eq. (ITBl) . This is 
done in Appendix O 

III. ENERGY OF CNOIDAL WAVES 

We next turn calculate the total energy of the Toda ring with a cnoidal wave. As a prelude, 

let us recapitulate the results of the trivial harmonic lattice, with the Hamiltonian in Eq. ([3]) 

truncated at the quadratic level. If we assume a phonon displacement Un = uq sin{kn — Ukt), 

then Uk = A/^|2sin||. The cycle average of the total kinetic energy and potential energy 

at wavevector k are easily found to be 
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KE = PE = A^M^Ksin^-, 



p = ^ = -NMul = -Nab\2sm-\^. (35) 

ul 2 2 2 

We next calculate the kinetic energy of the Toda ring using the cnoidal wave solution 
Eq. (!23|) as a Fourier series. Let us start with the kinetic energy expression: 

KE = ^j:(u,f 

j 

j 

00 2n L 

KE = 4 AT f ^ col) y , \ ,, sin^ n-. (36) 

n=l ^ ' 

Here we used Eq. ([7]), and the last line is obtained by taking the time average over a cycle, 
with the displacement from the series in Eq. ( 123|1 . 
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The potential energy of the Toda ring, averaged over a cycle, is calculated in Appendix 
[D| where we show that 



PE = A^- 
b 
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We expand Eq. ( ET^ , as in Eq. ( 1T^ in terms of the nome "q" to lowest order and find 



PE = N- 
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We now compute the ratio of kinetic to potential energies, and express it as a series in the 

nome: 

S = l+4g^in^(^) + 0(g^), (39) 

which is unity for small anharmonicity, as expected, and increases above unity for greater 
anharmonicity. 

In fact, we shall mainly compute the Fourier series for the mean square displacement 
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where ^{u) = J^^ ^Z{v)Z{u + v) [32]. Combining Eqs. ^ and Eq. (EE]), we find the total 
energy to mean square amplitude ratio 



KE 1 



P = 



{NMcoD 



n=l (l-q^n 



(l_q2n)2 



• 2 k 

sm n| 



w; 



TOO 1 q-^" ■ 2 k 

i;;^ (i_gSn)2 sm n^ 



(41) 



with Eq. ( IT6ll defining the spectrum U^. 

We can express the average kinetic energy in a more compact form using the expression 
in Eq. fITT]) for the displacements. Let us write 



KE 






Ui 



^^ r' dt {d,i<p - k) - d,me^''^*-'^-'^^^^^ 

Na 



2K 



du 
2b "'^ ./o 2^ 



( dn\ ^ ;^ ) ^°S 



1 + afc (dn^n 



AT' 



(42) 
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where we have set ak = ^1^ and used Eq. ( jDip . This expression is particularly useful 
if we want long wave length results since we can expand in a^ and find the exact elliptic 
parameter (i.e. m) dependent coefficients. The leading term for long wavelengths (up to and 
including 0{k^)) is 



KE 



Na 
'2b 

Na 



-«!- 



K 



du 



E' 



j,dn\ ^, 



26 



-ai 



+ 0{k' 



.E. 



2 E 1 ^ 

-(2 - m) (1 - m) - i—Y 



+ 0{k' 



Na 
'2b 



-ai 



~ < 



Na 2 



m 



m 



1024 



4 



1_6_ 
-m 



(m -)■ 0), 
(m -)■ 1). 



(43) 



31ogT^ 

=5 1— r 

We can also extract the long wave length behavior of Eq. (HOj) using [20'] as 



9 °° ^2n 

n=l ^ 



q' 



+ 0{k'), 



1 1, ,^2 ^;;^ 

- + - 2-m — + -^ 



1*^ 
1*^ 



24 6 
2 



256 
(log 



{l-m)f 



247r2 



(m -)■ 0), 



(m -)■ 11 



(44) 



We finally put together the results for the kinetic energy and the mean square amplitude to 
form the ratio 



p = -N{ab)k'^ 



256 



'm + m ) + 0(m 



p = N{ab)e^^^\og- 



{m -)■ 0), 
(m -)■ 1). 



(45a) 
(45b) 



27r2 " 1 — m' 

We see that the small m limit of Eq. (I45ap gives the long wavelength limit of the harmonic 
lattice result given in Eq. ( I35l) . In the opposite, strongly anharmonic, limit the kinetic energy 
grows logarithmically as ?7i — ;■ 1. 



IV. CONSERVED CURRENTS OF THE TODA RING 

We summarize in this section the construction of the conservation laws of the Toda lattice, 
and write out explicitly the first few conservation laws. For brevity we set m = a = 6 = 1. 
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The work of Henon [6] and Flaschka [TJ gives us a construction of the conservation laws 
starting with the Lax equation L = [L, A], where L and A are N x N matrices with entries 

^ _ „-rX _,p-|("m-"'i) _ A _,p-|("n-«m)l 

A„,^ = -V2 {5™,„+ie-5(«™-"") + ^^.^-le'^^""-"™)}, (46) 

so that the Lax equation leads to the original equation of motion, Eq. ([2]). Further, we can 
construct A^ independent conservation laws by taking the trace of the first A^ powers of the 
Lax matrix. With the definition 

J„ = ^TrL"; \<n<N (47) 

n! 

it has been shown that these are in mutual involution [l|, la, l7|, i.e. their Poisson brackets 
are all zero, ^Jn-,Jm\ = 0, and thus they form a complete set of independent conservation 
laws. We list the first few conservation laws: 

Ji = Ptotai, (48a) 

J2 = // + const. , (48b) 

J^ = lj2Pn + lj2P- {e(""-""+^) + e(""--«")} , (48c) 

n n 

n n n 

_| y^ ^(un-Un+2) _| Y^ ^2{un-u„+i) ^ (48d) 

n n 

We make extensive use of these conservation laws in later sections. 

V. PERSISTENCE OF ENERGY CURRENT AND ITS RELATIONSHIP TO 
CONSERVED CURRENTS 

In this section we study the decay and persistence of the energy current in the Toda 
lattice with periodic boundary conditions. The energy current is obtained from the energy 
density conservation law H{x, t) + dxJE{x, t) = 0, with a suitable discretization of the spatial 
derivative. To obtain the energy current, we write H = J2jP'j/i'^^) + i^ij^J' "W'here 
Vij = V{ui — Uj) = ^Q-Ku^-'^iHhj) Q^^(i r]{i,j) = {j — i)6i^\i-j\. The constant and linear term 
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in Eq. ([T]) can be omitted safely since they do not change the current. Thus the force on 
atom i due to atom j is 



Fij = -du.Vij = ar](i,j)e' 



-b{ui-Uj)r){i,j) 



We note the symmetry Vij = Vj^i and Fij = —Fji. Thus we write H = Y2- Hi with 
Hi = pf/{2M) + i J2j ^iji ciiid therefore denoting Vi = Pi/M, the rate of change of local 
energy is given by 

Hi = ViPi--Y^ Fij {vi - Vj) = 2 5Z ^*'J' *^^* + "^j)' 
j j 

= -Fi,,+i {vi + fi+i) - i^_i,i {vi + fi„i), 

= -JE{^ + l) + JE{^)■ (49) 

Therefore we may write alternate expressions for the energy current that are equivalent: 

Je = -^^Fi^i^iivi + Vi^i), or Je = -'^Vi{Fi_i^i + Fi^i+i), (50) 

i i 

with Fi-i^i = — ae"^*-"'"'"'"^''. We use the second form of the current above, and will set set 
M = a = b = lm the Hamiltonian Eq. ([T]) in the sequel. 
We compute the correlation function 

Cj^t) = {Mt)jEm, (51) 

where the average is taken in the canonical ensemble. This function does not decay to 
zero at long times but rather reaches a finite value. This is the phenomenon of temporal 
persistence, and is related to the integrability of the underlying model. The implication of 
this non vanishing of C{t -^ oo) is that the Fourier transform of C{t), namely the Kubo 
thermal conductivity, contains a Dirac delta function at zero frequency 6{uj). Thus the 
thermal conductivity has a Drude term in common with many other integrable models 
studied in recent years. 



Inspired by Mazur 2l| , we assume that the time integrated energy current can be written 



as a linear combination of all the conserved quantities of the model /„ 

ft 

t— >oo 



1 /"* 
lim - / JE{t') dt' = V aJn , (52) 
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for fixed coefficients an- Both sides of this equation depend on the initial state of the system. 
Again following Mazur, we argue that it holds, with the same a„, for almost all initial states 
on a constant energy surface. Here we will test Eq. f l52p for different initial states from a 
canonical distribution. The set of constants of motion /„ are not spelled out in detail in 
earlier work, although one expects that these are functions of the independent currents J„ 
in Eq. (jIT]). 

Using Eq. fl52|) . and a few chosen conserved currents, one can obtain a lower bound to the 



persistent part of the current current correlations 2l|] . The nice thing about this result is that 



even a single non trivial conserved current could help establish the existence of the temporal 



persistence ^2|. While in most applications of this idea, one has to be content with the 
result as a bound, one would also like to test the idea of completeness, i.e. to see if the bound 
is saturated in a case where one has a full knowledge of the conservation laws. A natural 
expectation is that if all the conserved currents that go into this expansion are known, then 
we should be able to get the exact value of the persistent part of the correlations. This is the 
operator analog of expanding vectors in a complete basis. Mazur gives the example of the 



XY model of magnetism in one dimension 231], where all the higher conservation laws are 



known in the quantum theory, since the problem is reducible to free fermions. In that case, 
he points out that exact persistent part of correlations can be obtained from a knowledge 
of these conservation laws. Initially, we will test the hypothesis that the set of conservation 
laws coincides with the constants of motion found in Eq. (l47l) . i.e. we begin by assuming 
that In = Jn- Later, we will discover that one needs to expand the set J„ to include further 
terms, such as bilinears in the J„. We begin by determining the a„ in terms of thermally 
averaged correlation functions. Multiplying Eq. (152|) (with /„ — > J„) by one of the conserved 
currents Jm and averaging over initial states with a Boltzmann distribution gives 

{JeJiti) = / ^ Q'n\JnJm) i (53) 

n 

where we used the fact that {JE{t)Jm) is independent of t. The a„ are therefore given by 

an = Y,{C-^)^,{JEJi), (54) 

I 

where C is the matrix of correlations of the conserved currents, 

Cnl = {JnJl) ■ (55) 
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Squaring Eq. (l52l) . and using Eqs. (153|) and ( 15^ . gives 



lim Cj^(t) = lim(Js(to)J£(to + t)) 



lini 

t—>ca 



tjo 



JEit') dt' 

n,m,l,k 

Y,{JEJi){C-'),,{JEJk), (56) 



k,l 



where we used ^^ Cmn C^i = (^m/ to obtain the last hne. Noting that Je is odd under time 
reversal, only the odd conserved currents, ra = 1, 3, 5, ■ ■ ■ , A^ — 1, in Eq. ( I48l) contribute. 

We have investigated the validity of Eq. fISB]) numerically for sizes A^ = 4, 6 and 8. 
All the numerical results in this paper are for parameter values M = a = h = 1 m the 
Toda Hamiltonian, Eq. ([1]). We prepare a large number (several thousand) of initial states 
appropriate to a temperature T (which we take to be T = 0.5) using standard Monte Carlo 



methods 



24l |. Starting from each of 



;hese states we perform molecular dynamics using 



a fourth order symplectic algorithm (25[, which combines high accuracy with long-term 
stability. We use a time step At equal to 0.05. 

Our results are presented in Figs, HI [5] and |6] for sizes A^ = 4, 6 and 8 respectively. The 
solid line is the data for Cj,^ (t) and the (blue) dotted line is the RHS of Eq. fl56l) including 
all the (odd) currents. While the long time limit of Cj^ (t) is close to the RHS of Eq. (!56|) 
there is a clear discrepancy. 

The discrepancy is removed if we note that the Lax currents J„ in Eq. (I48l) are not the 
only conserved currents but, in addition, products of these currents are conserved. This 
corresponds to enlarging the set of currents J„ in Eq. fl52|) to include bilinears JnJi- 

Thus the leading correction to Eq. (!56|) will involve products of two currents Jni = J„Ji 
where one of the indices must be odd and the other even for J„; to be odd under time reversal. 
Hence the total number of conserved "currents" to be included in the sums in Eq. fl56p is 
N/2 (the odd J„) plus (A^/2)^ the quadratic combinations with odd symmetry. Our results 
including the quadratic combinations of currents are shown by the (green) dashed line in 
Figs, m [5]and[6l The agreement is now excellent. We thus see that the set of currents that 
are involved in the expansion must include not just the Lax currents, but also products of 
these. We see that for the purpose of expanding an operator as in Eq. ( J52l) . the currents J„ 

20 



u 



0.36 - 
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0.34 - 



0.33 




FIG. 4: The correlation function of the energy current, defined in Eq. (I5ip . as a function of time 
for A^ = 4 particles. The (blue) short-dashed line is the RHS of Eq. ()56p including all the (odd) 
Lax currents. There is is a significant discrepancy with the long time limit of Cj^{t). The (green) 
long-dashed line is the RHS of Eq. (j56p including, in addition, all (odd) pairs of Lax currents. The 
agreement with the long time limit of Cj^{t) is now excellent. 



AJ| 



) 



(57) 



are not a linear basis, but rather Je seems to be algebraically dependent on the J„. Within 
the numerical scheme it is hard to determine if we have really saturated the persistent part 
by including just first and second powers of the Lax currents, but the remainder, if any, 
must be very small indeed. 

We also test Eq. fl5^ for the case of A^ = 4 by performing a least squares fit to minimize 

Nt 

Nt 

with respect to the a„, where tm = mAt, and Nt = t/At is the number of discrete times in 
the simulation. Including just the N/2 odd currents we find AJ|; = 0.016886, which is not 
extremely small and represents the difference between the dotted line and the long-time limit 
of the data in FigsHHHl However, including {N/2)'^ quadratic combinations of currents that 
project on to the energy current we get a much lower value, AJ"^ = 0.000220, as expected 
from the the good agreement between the dashed line and the long-time limit of the data 



1 ^* 
lim -— y^ JE{tm) - y^ a-nJu 



m=l 
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FIG. 5: The same as Fig. d but for A^ = 6. 
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FIG. 6: The same as Fig. [4] but for A^ = 8. 

in Figs HHHl As a consistency check, we compare in Table [T] the fit coefficients obtained by 
minimizing Eq. ( 157|) with those from Eq. ( 154|) which used equal time Monte Carlo results. 
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TABLE I: A comparison, for A^ = 4, between the coefficients a„ in the expansion of the energy 
current in terms of the conserved currents, see Eq. (|52p . The "time-series" results were obtained 
by minimizing AJ^ in Eq. (|57p . The "Monte Carlo" results were obtained from Eq. (j54p . 



time series 
Monte Carlo 



ai 



as 



0.4965 -0.9626 
0.4628 -0.9412 



«12 



ai4 



032 



034 



0.1262 -0.0591 0.0015 0.0178 
0.1358 -0.0681 -0.0042 0.0227 



The agreement is good. 

Mazur clarifies that the above expansion Eq. ( 152|) should be valid for almost all initial con- 
ditions on the surface of constant energy. Here we have seen that, in addition, the expansion 
works to high numerical accuracy in the more general case of a canonical distribution. 

VI. PERTURBED TODA RING AND DECAY RATES OF PERSISTENT CUR- 
RENTS 

We have also studied numerically the decay of persistent currents when a small perturba- 
tion away from integrability is added to the Toda Hamiltonian. The technique is the same 
as that described above in Sec. |Vl 

We consider two different perturbations: 

N 

3 



An 



'^ \^ t x3 



n=l 

N 



(58) 



y^ {Un - Un+lY 



n=l 



The temperature is taken to be T = 1, the time step to be At = 0.02, and the lattice has 
A^ = 64 particles with periodic boundary conditions. We have verified that reducing At or 
increasing A^ did not make a significant difference to the results. 

We compute the normalized correlation function of the conserved current J3, see Eq. (148 cl) . 
defined by 



Csit) 






(59) 



In the absence of any perturbation which breaks integrability, w or f in Eq. fl58|) . C^lt) 
is equal to unity. It is also equal to unity at t = for any Hamiltonian. Figure [7] shows 
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V(x) = exp(x)-x-l + wx^/3 




FIG. 7: The normalized correlation function, defined in Eq. ()59p . for the current J3, for several 
values of the cubic perturbation w in Eq. (I58p added to the integrable Toda Hamiltonian. The 
size is A^ = 64. Including the perturbation, the system is no longer integrable, so the correlation 
function tends to zero on a timescale which diverges as w ^ 0. 



data for C^it) for the cubic perturbation in Eq. fl58|) . for several values of the strength of 
the perturbation w. As expected the J3 correlation function decays to zero on a timescale 
which decreases with increasing w. 

We assume that the data fits the scaling form 



Cs{t) = Csitw'^) 



(60) 



where is a crossover exponent indicating that w is a "relevant" perturbation for the 
integrable model. Figure [H] shows scaled data for the cubic perturbation in Eq. (!58l) assuming 
= 1.9 which gives the best data collapse for the largest sizes. 

However, as shown in Fig. [9] it is also possible that the data for even larger sizes will 
collapse for = 2. The result = 2 (in a quantum version) follows from Fermi's golden 
rule. It is commonly encountered in quantum integrable systems 22|, where the current 
decays due to the addition of a term V in the Hamiltonian that destroy integrability. The 
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FIG. 8: A scaling plot, according to Eq. (pO]l . of the data in Fig. (Tj with crossover exponent (p = 1-9 
which gives the best data collapse for small w. 

Fermi golden rule states that the decay rate of a state i is given, to lowest order, by 

1 9~- 

7 = tEk^i^i^')I'^(^^--^0- (61) 

This may be interpreted as the decay rate of a current or a quasiparticle that is infinitely 
long lived in the integrable case. Hence we expect that an integrable current analogous to 
J„, and its correlations would decay in this fashion, whereby we expect = 2. 

We note that the data in Fig. [9] only collapses for very small values of w and very large 
times, indicating that "corrections to scaling" are large. 

We have also studied the effect of the quartic perturbation in Eq. (!58|) and show the 
unsealed data in Fig. [TOl We expect the same crossover exponent (p for the quartic pertur- 
bation as for the cubic perturbation, and indeed the data is consistent with the asymptotic 
exponent being = 2 as shown in Fig. [TTJ 
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FIG. 9: Same as Fig. [8] but for (p = 2, which is the expected result as discussed in the text. It is 
plausible that this value of will give good data collapse in the scaling limit ii; — )■ 0, t — J- oo. This 
asymptotic regime is reached only for very small values of w (and correspondingly very long times) 
indicating that corrections to scaling are large. 

VII. CONCLUSIONS 

In this paper we have shown that the result of Toda for the frequency of the cnoidal 
waves of the Toda lattice is only correct for weak anharmonicity or long wavelength. We 
give a general expression for the frequency, as well as results for the average kinetic energy 
and mean square displacement in the Toda ring. The distinction between the dispersion 
relation of Toda, Eq. flB7|) . and the one found here, Eq. flT6l) . is important only in the limit 
of large wave vectors and high anharmonicity-in all other cases it is negligible. 

In addition, we have discussed the conserved currents of the Toda model in some detail. 
In particular, we showed numerically that the persistent part of general currents can be 
expressed in terms of the conserved currents according to Eq. flS^ . provided one includes 
not only the Lax currents, but also quadratic combinations of the Lax currents (which are, of 
course, also conserved). Finally, we have studied the decay of the conserved currents when a 
perturbation is added to the model which destroys integrability. The timescale for decay is 
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FIG. 10: The same as Fig. [7] but for the quartic perturbation v in Eq. ([58 
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FIG. 11: A scahng plot with <p = 2, as in Fig. [Dj but for the quartic perturbation v in Eq. (j58|) . 

governed by a crossover exponent 0. Our numerical data is consistent with the value = 2, 
which can also be obtained by Fermi golden rule-type arguments. 
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Appendix A: Elliptic Theta Functions and Integrals 

The most convenient version of these functions is given in Whittaker and Watson, Chapter 



XXI 



18| and Abramowitz and Stegun 26| . We note here the definitions of the required theta 



functions and their Fourier series: 

oo 

q ^ ^-^K'/K^ Q^^^^ = 1 + 2 ^ (-1)" g"' cos(2n;2), (Al) 



n=l 

oo 



K = K{m) = dx =, ei{z) = 2 J2 (-1)"?^"+^^' sin(2n + l)z 
Jo \/l — m sin^ x n 



V 1 — m sm" X 
K' = K{1 - m), e2{z) = 2 J2 ^^"^^^' cos(2n + l)z, 



n=0 

oo 



n=0 

oo 



E= I dx Vl-msin^x, ^3(z) = 1 + 2 V g"' cos(2nz). (A2) 

•^0 n=l 

We note that for small anharmonicity, all the elliptic functions can be expanded equivalently 

in terms of the parameter m or the nome q. These expansions are related since we may 

expand 

q= — m -\ m -\ m -\ m +(J{m). A3 

^ 16 32 1024 2048 ^ ^ ^ ^ 

We note that as m — > 1, K{m) ~ log{~7T^=} and E{m) ~ ^. Further, a.sm — )■ 1, the theta 

functions are more easily calculated by using the Jacobi transformation [18i] (p 475), where 

we display the parameter m 

9n{z\m) = c„ -^e"^^'^^ Gn{ij^z\m), (A4) 

with ci = —i, and the remaining Cj = 1. For computing the soliton spectrum, we will need 
the expression for n = 1 

^i(z|m) = 2 — e-^"'^> ^(-l)"e-"#^("+^)'sinh(2n + l)—z. (A5) 

n=0 
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Appendix B: Toda's cnoidal wave frequency calculation, duality and boundary con- 
ditions 

Consider N atoms in 1-d, with displacements Ui,U2. ■ .uj\f, which we call the bulk dis- 
placements. Additionally, there are two "boundary" atoms, atom to the left of atom 1, and 
atom A^ + 1 to the right of atom N. Here we do not assume periodic boundary conditions 
but instead consider two different types of boundary conditions: 

1. "clamped" boundary conditions for which uq = un+i = 0, 

2. and "open" boundary conditions for which, it turns out, we will need uq = "Ui and 

UN+1 = UN. 

We will actually focus on "mixed" boundary conditions, i.e. free at one end and clamped at 
the other. 

The interaction between two atoms is represented by a nearest neighbor term V{uj — Uj-i) 
with a suitable function V, the exponential interaction as in Eq. ([T]), or a harmonic term for 
comparison. The Lagrangian is 

, N Jr 

i=i j=Ji 

where Ji = 1 (2) for clamped (open) boundary conditions at the left, and J,. = A^ + 1 (A^) 

for a clamped (open) boundary conditions at the right. Even for open boundary conditions, 

the boundary sites do not have any kinetic energy since they are "fictitious" , i.e. are simply 

there to impose the necessary boundary conditions. 

We introduce the particle separations Tj = Uj — Uj^i as new generalized coordinates, these 

and their canonically conjugate momenta Sj found below, are the "dual variables" of Toda. 

It follows that the inverse relations are Uj = X]!?=i'"'- '^^^ kinetic energy is expressible 

in terms of the time derivatives, Tj. The canonically conjugate "momenta" to the Tj are 

Sj = ^ = M '^i^j ui, whereby, for 1 < j < A^ — 1, we get Miij = Sj — Sj+i, and sn = Mun- 

Thus we obtain the Hamiltonian in terms of the dual variables as 

^ = ^E(^.-^.>i)' + ^4 + E^(^.)- (B2) 

i=i j=Ji 

Hamilton's equations of motion are therefore now 

1 dVir) 
f, = —{2s, - s,_i - s,+i), s, = ^ for 2 < J < AT - 1. (B3) 
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The case of a clamped right boundary is comphcated to deal with so we will always choose 
it to be open. We therefore have Jr = N, and the kinetic energy term of the A^*'^ atom is 
written as ^{s^ — stv+i)^ with Sjv+i = as the boundary condition at all times, so that 
the equation Eq. ( 1B3P is extended to j = N. 

On the left side boundary, the clamped case with Ji = 1 can be dealt with easily by 
extending Eq. ( 1B3|) to j = 1 by introducing an sq and requiring that Sq = Si so Eq. f IBSP 
has the same form for j = 1 as for j > 1. The open case on the left boundary has a missing 
V{ri) in the potential energy so that ri is a cyclic coordinate (so fi = 0) and we must 
require Si = constant. 

For the case of the harmonic chain, these boundary conditions are easily imposed on 
the solution e^*'^** cos{krj + S) and one reproduces the various integer and half integer 
quantization of the wave vector k. We note that the solutions are not traveling waves, 
but rather products of functions of space and time, since in all cases here we must set the 
boundary variables to be time independent constants. 

For the Toda lattice, we choose V{rj) = | (e"^^'^ — 1 + bvj) so that the equations of 
motion Eq. ( JBSp are written for the typical case of left-clamped and right-open boundaries 

s ■ b 

l_ . = jT [sj-i + Sj+i - 2sj] , 1 <Sj < N, and Sq = Si, sn+i = 0. (B4) 

The remarkable insight of Toda in solving this equation, was apparently inspired by seeing 
an addition identity of the Jacobean zeta function Z[u] that we met earlier in Eq. fl20l) . 
The salient features of this function are as follows: periodicity Z[u + 2K] = Zju], parity 



Z[—u] = —Z[u], and hence nodes at Z[0] = = Z[K]. The relevant identity is 15| 



Z[u + v] + Z[u -v]- 2Z[u] = f ["] , ■ (B5) 

Comparing Eq. ( 1B4I) and Eq. ( JB5I) . one sees that these are very similar, provided we make 
the hypothesis that Eq. ([5^ should collapse to an ordinary differential- difference equation. 
Thus we are obliged to combine the space and time dependence into a single variable 0j = 
kj — ut + S, and this helps in solving the Eq. (]B4p for the bulk. We thus can only solve 
for traveling waves. However, it makes it impossible to satisfy the boundary conditions, 
since the latter involve time independent vanishing of certain constants. With Toda, we will 
ignore the boundary terms and write down the solution for Sj that is implied by Eq. ( JBSp . 
With a scale factor -^ we relate u and as u = ^ K, so that increasing u by its natural 
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periodicity 2K winds the phase by 27r. Then we see that ^ = —^^-^^ and j — )■ j + 1 
increases u -^ u + — whence v = —. The mapping is complete with a scale factor relating 
Sj to Z[u], and determining a; as a function of k gives the two Toda solutions: 

.,(t) = Ty (^) Z[^ikjTult + 6)] (B6) 

T ah Tx 1 ,^„, 



Using Miij = Sj — Sj+i and the expression Eq. fl2UI) . and integrating once we find the 
displacement in terms of the theta functions: 



_1 e,{likjTu;lt + 5)) 



^. = T log . .wV. T^ ' r. ^^^ . (B8) 



in agreement with our Eq. f fTOj) . 

In Sec. IllCl we showed that the same functional form as Eq. (]B8|) describes oscillatory so- 
lutions with periodic boundary conditions, but with a different expression for the frequency, 
Eq. (IT6l) . in place of Eq. (IB 6 1) . This situation requires a few clarifying remarks. 



Eq. ( ]B6p is in a traveling wave form, and since we cannot superpose two non linear 
waves, e.g. with the two signs of the time dependence, this prevents us from satisfying 
the various boundary conditions (clamped and open) discussed above. Thus the Toda 
dispersion relation is not the solution of the problem he starts with. 

The spectra in Eqs. ( IT8|) and (TT6|) are very close for small wave vectors, or for small 
values of the parameter "m", i.e. for weak anharmonicity. Fig([T]) illustrates the two 
dispersions. 

A reader might wonder if Eqs. (IB7I) and flT6|) are not actually identical, with the help 
of some obscure identity, e.g. in Ref. [l8|. However, this cannot be the case for two 
reasons. One is mathematical. While the expression in Eq. (IB6P is doubly periodic in 
the complex k plane, the expression in Eq. ( IT6l) is not-the 6 functions have periodicity 
factors attached for translations along the imaginary axis. The other reason is physical. 
If the two expressions were somehow identical, the potential energy average Eq. (ID3P 
would be identically zero for the periodic system, which is impossible since V^(x) > 
with the equality only at a; = 0. 
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Toda is correct to take the functional form in Eq. f lB6P seriously. However for periodic 
boundary conditions, the correct frequency of the cnoidal wave is not the one in Toda's 
works ij, but rather Eq. f ll6p . which appears here for the first time, to the best of our 
knowledge. Toda's solution is only correct in the limit of weak anharmonicity or long 
wavelength. 



Appendix C: Alternate derivation of the solution with periodic boundary condi- 
tions. 

In this appendix we give an alternative derivation of the dispersion relation of the Toda 
ring in Eq. f lT6|) by substituting the Fourier series in the second line of Eq. f l23|) into Eq. (|9]) . 
Taking the second derivative of this with respect to we obtain 



oo 



LHS Eq. © = 4 a;2 ^n -^ sinn(0 - -) sinn-. (CI) 

n=l ^ 

To obtain the expansion of the RHS of Eq. (jH]) we start with Eq. fl25]) and use the Fourier 
expansion of sn^ due to Jacobi [19] (again!) 

sn ( x) = —{l-—} ^V^i ^ cos2na;. (C2) 



1 - ^ 

ra=l 



Thus we find 



RHS Eq. B = ^=^ ^ g „ ^ s,„ „# - -) s,n „-. (C3) 

The series in Eqs. (ICip and (]C3|) are seen to be identical, with the choice of the dispersion 
in Eq. (|T6|). 



Appendix D: Potential energy of the Toda ring 

Let us rewrite the potential energy term Eq. (fTTl) . with representing any one phase 
factor; 
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^bu„ — bun+i 



,dk{<p)-dk{<t>+k) 



1 em 






4V2. 






u 



,K^ 



k ^2 ] ^^2(kK^ 



m sn 



,i^ 



TT 



(Dl) 



We used the two dispersion relations for the chain and the ring, Eq. flT6l) . and Eq. f lT7|) to 
proceed in this equation. 

We now recall = ujkt and average this expression over a single cycle in time, (i.e. 
L'" ■ ■ ■ jr "where T^ = ^), or m = -^0, with < m < 2/^. We ignore the site index, since 
each atom has the same average value. We use 



2K 



1 fE 



—- / dusn {u) = — ( — -I 
2K Jo m \K 

so that the time average of this term can be written using the expression Eq. flB7p as: 

— \ 2 



gdfc(fli)-dfc(0+fc) ^ :^ 






and thus the potential energy average over a cycle is 



PE = A^- 





— \ 2 

col 



(D2) 



(D3) 
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TTU vr ^ sin^ 






2KK' 2K ^ sinh^' 
Since this is in the form of a Fourier "sine-series" , we use the "sine-transforms" version of the 
Poisson summation formula: 

oo oo „oo 

^ 5r(n)sin — = uq ^ g{u- 2Truou), where g{x) = / dy sm{xy) g{uoy), 

n=-oo ^0 j.=-oo -^-"^ 

with uq = K/n and g{n) = ^ sinh(nK'n/K) ' yielding the required result. Similarly Eq. (j23|) 
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